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^ ! ABSTRACT 



04 ' Context. Determining the evolutionary stage of a Young Stellar Object (YSO) is of fundamental importance to test star formation 

' theories. Classification schemes for YSOs are based on evaluating the degree of dissipation of the surrounding envelope, whose main 

effects are the extinction of the optical radiation from the central YSO and re-emission in the far-infrared to millimeter part of the 
electromagnetic spectrum. Since extinction is a property of column density along the line of sight, the presence of a protoplanetary 
disk may lead to a misclassification of pre-main sequence stars with disks when viewed edge-on. 

Aims. We performed radiative transfer calculations to show the effects of different geometries on the main indicators of YSO evolu- 
^ , tionary stage. In particular we tested not only the effects on the infrared colors, like the slope a of the flux between 2.2 and 24 fjm, 

^ ■ but also on other popular indicators of YSO evolutionary stage, such as the bolometric temperature and the optical depth of silicates 

^ ' and ices. 

' Methods. We used the axisymmetric 3D radiative transfer codes RADMC and RADICAL to calculate the spectral energy distribution 

including silicates and ice features in a grid of models covering the range of physical properties typical of embedded and pre-main 
sequence sources. 

04 ' Results. Our set of models compares well with existing observations, supporting the assumed density parametrization and the adopted 

^ ' dust opacities. We show that for systems viewed at intermediate angles (25° - 70°) the "classical" indicators of evolution are able to 

0^ , classify the degree of evolution of young stellar objects since they accurately trace the envelope column density, and they all agree 

CO . with each other. On the other hand, edge-on system are misclassified for inclinations larger than ~ 65° ±5° , where the spread is mostly 

due to the range of mass and the flaring degree of the disk. In particular, silicate emission, typical of pre-main sequence stars with 
disks, turns into silicate absorption when the disk column density along the line of sight reaches 1 • IO^"cm"-, corresponding e.g. to 
a 5 ■ 10"^ Mq flaring disk viewed at 64°. A similar effect is noticed in all the other classification indicators studied: a, Tboi, and the 
H2O and CO2 ice absorption strengths. This misclassification has a large impact on the nature of the flat-spectrum sources (a ^ 0), 
whose number can be explained by simple geometrical arguments without invoking evolution. A reliable classification scheme using 
' a minimal number of observations is constituted by observations of the millimeter flux with both a single dish and an interferometer. 
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, 1. Introduction extinguished and scattered by the suiTounding envelope. It is 

■ - - ■ evident that measuring the "relative intensity" of the two com- 

A correct classification of the evolutionary stages of young stel- ^^^^^^^ ^^uld furnish the contribution of the envelope to the 

lar objects (YSOs) is fundamental to understand the time-scales ^j^^j^ ^y^^^j^ j^^^g^^ low-mass protostellar systems were orig- 

involved in the star formation process. Since protostars grad- j„^jjy classified according to the amount of extinction towards 

ually accrete and disperse their envelope material during their j^e protostar, using the ratio between near (2 ^ m) and mid- 

hfetime, the amount of residual envelope material can be used to infrar ed fluxes (25-60 ^m ) (a= dlo2(AF ,)/dlo 2(A). \ Adams al\ 

measure the degree of evolution of the protostellar object. The |T987l:lGreene etd][lW4 ). or by the relative amount of emission 

presence of a surrounding envelope has several characteristic ef- f^^j^ j^e envelope, evaluated by the average frequency of emis- 

fects on the observables of the YSO which can then be linked • ., r ir j , r ip j \ »j» uiiuj* 

, . ^ , T • , . sion «v) - vFydv Fydv) converted to a blackbody temper- 
to the evolutionary stage or the system. In particular, the enve- . "rj c t jJlinmK t-u u- / u 

rr , , ,. ., . .or-.^^ r aturc (T\„\, M vers & Laddlll993h . The combination of observa- 

lope strongly affects the spectral energy distribution (SED) of ^^^^^ ^ three-class classification with class 

the protostellar system, extinguishing and reprocessing the pro- ■, » ,.u 1 ^ ^ a » » ,.u ■ »f 

^„ ,. . ™ . , „ . , , ° . . I representing the least evolved protostars with a significant frac- 

tostellar radiation. Thus, instead of a single blackbody emission, ^- e »n 1 ■ 1 a u ■ a 

^ ' ° , . ' tion of the mass still left in the envelope and showing a>\) and 

the SED of a youiig protostellar system appears as a combina- ^,^^^^^5^ j^. ^j^^^ corresponding to the classical T Tauri stars 

tion of two thermal components: one corresponding to the enve- -^u a- \ a • ■ e a ■ ■ j o , ^ a 

. . ,Xnrr^ , , , , . , , r with disks dominating the infrared emission and with -2 <a< 

lope emission (J „ < 100 K) and the other that is the result of K<rboi<2800 K; and class III pre-main sequence stars 

the emission of the central protostellar object (T^g > 2000 K) ^^.^^ ^^^^ p^^^^^^^ circumstellar disks only in the mid- 

infrared and have a< -2 and rboi>2800 K. The earlier class 
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deeply embedded protostars that are invisible in the infrared. 
Finally, some authors refer to the sources with -0.3 <a< 0.3 
as flat-spectrum sources a nd view them as transitional objects 
between class I and II (e.g. jGreene et alJfT994l) . 

This classification has had an enormous success and a huge 
impact on star formation studies, but it is subject to misclassifi- 
cation. In particular, it does not take into account that the dissi- 
pation of the envelope is not a spherical process; in fact, while 
part of the material is dispersed in the outflow, the bulk of it is 
accreted into a disk, with the result that systems viewed edge- 
on would present extinction values typical of younger systems 
with more massive envelopes, thus leading to a misclassifica- 
tion. This effect is clea rly sh own in the theoretical w ork of, e.g., 
IWhitnev et all (l2003allH) and lRobitaille et alj (l2006h . To empha- 
size the difference between the observational classification and 
the intrinsic properties of the YSO, Robitaille et al. (2006) intro- 
duced the alternative nomenclature "stage I, II, III", equivalent 
to the Lada classes, but referring to the true nature of the source 
regardless of its observed properties. These authors distinguish 
stage I from stage II YSOs according to their accretion rate of 
the envelope, with the dividing line set at Menv = 10"'' Mq yr"' 
(stage II and III are defined according to the accretion rate of 
the disk instead of that of the envelope). In the equations for a 
rotationally flattened envelope, it easily can be seen that Menv is 
ultimately a measure of the volume density of the envelope at 
the centrifugal radius; for a 1 M© star with a 300 AU centrifugal 
radius Menv - 10"^ Mq yr ' is equivalent to an H2 number den- 
sity of 1.8- 10^ cm 

A possibility to observationally disentangle the extinction com- 
ing from envelopes to that produced in the disks can come from 
their very different temperature structures, with disks going from 
30 to 1000 K and envelopes from 10 to 200 K. With such dif- 
ferent conditions, dust grains are expected to be very different. 
Specifically, the envelope grains are expected to be coated by 
ice unlike the fraction of disk grains with temperatures above 
100 K. Thus, comparing the evolutionary stage inferred by the 
infrared colors, like a, with Tboi or the silicates and ice opti- 
cal depths could provide a tool to distinguish between stage I 
YSO and edge-on disks. Including the presence of ices (H2O 
and CO2) in our opacity tables enables us to predict not only 
the i nfrared colors, as effi ciently presented in the recent litera- 
ture (iRobitaille et al.ll2b06h . but also the optical depth of spectral 
features, in an attempt to broaden the number of observables for 
YSO evolution diagnostics. 

In this paper, we present model predictions for a series of 
YSOs covering a large range of disk mass, envelope mass and 
stellar luminosity, and exploring the effect of inclination on the 
emerging SED and spectral features. Section |2] presents the ra- 
diative transfer model we used and the adopted ice-coated dust 
properties; Section[3]compares the results of the calculation with 
existing data; Section|4] summarizes our conclusions. 

2. The radiative transfer model 

2.1. Physical Structure 

We used the three -dimensional axisymmetric r adiative trans- 
fer code, RADMC jDuUemond & Dominikl2004 . to predict the 
continuum emission of a grid of 63 low-mass protostellar sys- 
tems, each viewed at 7 different inclinations. The code takes as 
input the characteristics of the star, the density structure of the 
surrounding material and the dust optical properties, and out- 
puts the temperature structure a nd the scattering source function 
using a Monte Carlo technique (iLucvll 19991 ; iBiorkman & WoodI 



|200T1) . In this version of the code, scattering is supposed to 
be isotropic. Ray-tracing is then performed using RADICAL 
(Dullemond & Turolla 2000) to produce SEDs and images. 

Given the high density reached in the center of our models 
the exact spectral energy distribution of the star does not matter, 
because the emerging spectrum will be fully reprocessed. For 
this reason, we assume the (proto)-star SED to be a blackbody 
at 2500 K and we vary only its total luminosity in our grid be- 
tween 2, 6, and 18 Lq. 

The density structure is the sum of three different components: 
the disk, the envelope and the outflow cone. The adopted density 
structure for the disk has a power-law dependence along the ra- 
dial coordinate and a Gaussian dependence in height, and can be 
expressed as 

Xo(r//;o)-' f 1 [ rcosfl ]'] 

Pdisk(r, 0) = — — exp ^ - - > , 

V2^//(r) { 2 [ H{r) \ j 

where 6 is the angle from the axis of symmetry. The variation 
of scale-height, i.e the flaring of the disk, is described in the 
function Hi/) - r ■ HqIRq ■ (r/Rp)^^'', corresponding to the self- 
irradiated passive disk of IChiang & GoldreichI (Il997h . In our 
grid of models we fix Rq to be 200 AU, Hq equal to 40 AU, 
and we vary the total disk mass between 5 ■ 10"^, 5 ■ 10"^, and 
5 ■ 10-* Mo. 

The envelope density follows the theoreti cal structure fo r a rotat- 
ing and collapsing spheroid as derived bv lUlrichl(ll976l) . defined 
by 

Penv(r, 0) = PO 1 + -\ \- - + COS^0o , 

\ r I \ cos 6qI \2 cos Op r j 

where 6q is the solution of the parabolic motion of an infalling 
particle given by r / R^oticos Bp - cos 0)/(cos 00 sin^ 0o) = 1, ^lot 
is the centrifugal radius of the envelope, and po is the density in 
the equatorial plane at the centrifugal radius. The outer radius 
of the envelope is fixed to 10000 AU and the centrifugal radius 
set to 300 AU, while we vary po to cover the range of masses 
between 0.01 and 4.0 Mq in seven steps. The inner radius of 
the envelope is determined by the outflow cavity (see below), 
which crosses the equatorial plane at ~20 AU. Comparing our 
parametrization of the density structure to that of lRobitaille et al.l 
(j2OO60 and considering that our models have a fixed size and 
centrifugal radius, we can convert their assumed boundary for 
stage I/II of Menv - 10~^ Mq yr"' to a more intuitive bound- 
ary in envelope mass of Menv= 0.07 Mq. We note that this value 
would vary by factors of a few if different star mass, centrifugal 
radius and envelope outer radius are assumed. Given this spread, 
and considering that circumstellar disks more massive than one 
tenth of a solar mass would be gravitationally unstable, we chose 
to round the boundary value between stage I and II YSO to 
Menv= 0.1 Mq. This boundary envelope coiTesponds to column 
densities of molecular hydrogen between 1.8 and 6-10^' cm"^ at 
inclinations between 25° and 85°. 

We include a streamline outflow cavity by setting the density of 
the region where cos 0q is larger than cos 15° to the same value 
of the density of the envelope at the outer radius. This results 
in a funnel-shaped cavity which is conical only at large scales, 
where it presents a semi-aperture of 15°. The total column den- 
sity inside the outflow cavity can be found in Table [T] 

The output temperature structure and scattering source func- 
tion are then ray-traced at 7 different inclinations between 25° 
and 85° and spaced by 10°. The explored parameter space in 
Lstm, Mdisk, Menv, and inclination results in 441 models. 
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In Fig. [T] we show the temperature and density structure for 
two models with different envelope mass. The two structures 
are plotted on two different scales in the left and in the right 
panel to emphasize the disk and the envelope. The temperature 
of the most massive envelope is in general lower than that of the 
less massive ones. This effect is produced by the higher density 
present in the inner wall of the envelope resulting in a larger frac- 
tion of energy that is absorbed closer to the star. Note also that 
most mass in the disk is at temperatures lower than 90 K and thus 
its dust will be covered by H2O ice. For the same two models, 
in Fig. |2l we show the contribution of the disk and of the enve- 
lope to the column density along the line of sight as a function 
of the viewing angle. In the left panels the total column density 
is presented while in the central (right) panels the contribution 
of iced (bare) grains is isolated. The disk material dominates the 
column density for the lines of sight that are inclined more than 
60°-70° from the axis of symmetry. This is the case not only for 
the total column density but also for the column of iced grains 
in a pencil-beam . In fact for extreme edge-on systems the disk 
contribution to the column density of ices can be a few times 
larger than that of a massive envelope. 



2.2. The dust opacity model 

The optical properties of the dust grains used in this paper were 
obtained by averaging the properties of a distri bution of silicates 
and grap hite grains covered by ice mantles. Pontoppidan et al.l 
din prep.h obtained a large grid of opacity tables for silicate and 
graphite grains with different size distributions and different 
amounts of ice mantle and compared those to the extinction 
observations between 1 fim and 1 mm obtain ed by the "Cores 
to Disks" (c2d) Spitzer legacy program ( Evans et al.l l2003h 
and presented by Huard et al. (in preparation). The best fit 
to the extinction law of dense cores is obtained by a mixture 
of 71% silicates with a siz e distribution similar to that of 
IWeingartner & Draind (1200 ll) with an awD-paramete^3 equal to 
-2 and turnover point for radii larger than 0.3 /im covered by 
a mantle of ices with a water ice abundance of 3-10"^ relative 
to H2, and 29% of carbonaceous grains w i th a s ize distribution 
similar to that of IWeingartner & Draind (1200 Ih characterized 
by a turnover point at radius of 4.5 fim an awD-parameter of 
-2 and a water ice abundance of 1.510 relative to H2. This 
dust mixture was used in the whole grid with the exception of 
the regions with temperatures higher than 90 K where all the 
ices were removed. This was performed by running the Monte 
Carlo code once using only the opacities with ices in the whole 
grid to calculate the temperature structure, and then using this 
information to substitute the dust opacities in the grid cells that 
are found to have temperatures higher than 90 K; a new Monte 
Carlo run was then performed to determine the correct scattering 
function. Given the minimal differences in the opacities with 
or without ices the temperature structure is almost unaltered. 
Pure CO2 ice desorbs at a slightly different temperature than 
H2O (-65 K vs. -90 K in space, ISandford & Allamandolal 
[1991 : since we kept the relative abundance of the two ices fixed 
to 20%, the differential desorption in the temperature range 
65-90 K is not taken into account in the current grid. 



' The ttwD of IWeingartner & Draind ( 1200 ih (their notation) is a pa- 
rameter of the grain size distribution, not to be confused with the mid-IR 
slope of the SEDs 



3. Results 

3.1. Spectral energy distributions 

In Fig. [3] we present the spectral energy distributions of a se- 
ries of models with Mdisk=5 ■ 10"^ M0 and Lstar=6 Lq. Fluxes 
were calculated integrating the emission inside circular apertures 
with radii increasing from 280 AU for near infrared wavelenghts 
to 10000 AU for millimeter wavelengths, to simulate the cur- 
rent observing capabilities for sources at 140 pc, the distance of 
nearby star-forming regions such as the Taurus and Ophiucus 
molecular clouds. In the six panels we show six different incli- 
nations while the sequence of SEDs shows the variation of en- 
velope mass increasing from 0.01 to 4 Mq from top to bottom. 
The SEDs are normalized to be equally spaced at 8 fim in or- 
der to show more clearly the differences in the silicate features. 
Systems closer to edge-on show a higher degree of numerical 
noise because the determination of the temperature and of the 
scattering source function in some of the cells in the mid-plane 
is made difficult by the high densities that limit the number of 
photons penetrating the innermost cells. Systems closer to face- 
on, depicted in the top three panels, show the standard sequence 
of SED from stage I to stage II (from bottom to top), with the sili- 
cate feature turning from absorption to emission and the ice fea- 
tures disappearing. The turning point from absorption to emis- 
sion happens when the column density of molecular hydrogen 
drops below ~ 1 ■ 10^^ cm"^ in a pencil beam, corresponding to 
an envelope mass of at most 0.20 Mq in our models. The "evolu- 
tion" of the SED is measurable also from the change in a (calcu- 
lated from 2.2 fim to 24 fim, and reported in Fig. [3), turning from 
positive to negative with decreasing envelope mass, and from the 
bolometric temperature, increasing from 30 K to 1200 K. These 
results show, once again, that the Lada classification is accurate 
for face-on systems. More details on the models shown in Fig. [3] 
such as column densities, a, and Tboi, are reported in Table [1] 
In the three lower panels of Fig. |3] we start to see the effect 
of the disk in absorption with the disappearance of the silicate 
emission even in the bare disks and the corresponding negative a 
values. This effect is visible not only for the extreme 85° inclina- 
tion but also for the more moderate value of 65°. A finer grid of 
inclinations, performed on a model with a very small envelope 
(10""* M0), shows that the turning point for the silicate feature 
occurs at 64° for the 5-10"^ Mq disk, and at 73° and 58° for the 
5 ■ 10"^ Mq and the 5 ■ 10"^ M0 disks, respectively. This varia- 
tion reflects the different inclinations where a column density of 
~ 1 ■ 10^^ cm"^ is encountered in our models. We showed above 
that at the same column density a also turns from positive to neg- 
ative. This means that, even if we consider only the disks with 
the smallest masses, all the stage II sources seen in the solid an- 
gle between 70° and 90° , equal to 47r(cos(70) - cos(90)) or 34% 
of the total, may be mis-classified as class I or flat-spectrum 
sources. Recent large-scale mapping of the Perseus, Chamaleon 
II, and Serpens star for r tiing regions perfor r ned by the c2d tearn 
( lJ0rgensenet al.ll2006t iPorras et al.lf2007t iHarvev et alj l2007l) 
show that the ratio between flat-spectrum sources and class II 
objects ranges between 19% and 47% with an overall average 
of 26%. Our simulations can explain the observed number of 
flat-spectrum sources by a geometrical argument alone, without 
necessarily invoking evolutionary effects. 

3.2. Testing the classification schemes 

To better quantify the effect of inclination in the classification of 
YSOs, we plot in Fig.|4]four evolutionary indicators as a func- 
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tion of the mass of the envelope in our models: a, Tboi, the 10 fim 
silicate feature strength, and the 3 jim water ice stretch feature. 
The edge-on models (inclination equal to 65°, 75° and 85°) are 
marked with empty triangles offset by 15% in mass for clarity 
and the face-on models (inclination equal to 25°, 35°, 45° and 
55°) with filled circles. The colour coding separates models with 
different envelope mass and is thus a measure of the evolution 
from stage I to stage II YSOs. In all four panels, it can be seen 
that the face-on models form a well-defined sequence that allow 
us to correlate the envelope mass to the evolutionary indicators. 
In fact, no filled circle crosses the empirical boundary for class 
I/II for Man <0.1 in any of the four panels. This is not the case 
for the edge-on models, which appear less sensitive to the de- 
crease of envelope mass, resulting in similar values of a, Tboi, 
and silicate and ice features for both stage I and stage II ob- 
jects. As a result, edge-on disks can be erroneously classified as 
class I YSOs when using these evolutionary indicators. This is 
because a, T\,oi and the silicate feature are sensitive to the col- 
umn density along the line of sight, and the additional column 
of material present in the disk has a similar effect on the SED 
as a more massive envelope. However, even though the column 
density through an edge-on disk can exceed the column density 
of the most massive envelope by one or two orders of magnitude 
(see Table[T]|, none of the edge-on disks (Me„v<0. 1) present val- 
ues of a as high as those observed in the models with Menv= 
4 Mo. This is due to the fact that a large part of the 2 fj.m flux is 
scattered back into the line of sight through the material above 
the disk and this scattered light is never absorbed again, lead- 
ing to large 2 fj.m flux and consequently much smaller a values 
than those measured in massive envelopes, where the extinction 
at 2 fim comes from the sum of absorption and scattering. In 
other words scattering allows lines of sight with low absorption 
at 2 fim for edge-on disks unavailable for massive envelopes. 
The effect of inclination on Tboi (Fig.|3)) grows in strength from 
stage I to stage II YSO. We note that even for the least massive 
disk that we studied, at large inclinations Tboi crosses the limit 
for class I sources (650 K) and would lead to a classification 
as class I. Given the much larger number of observed class II 
than class I sources i n star-forming regions (392 vs 104 in the 
combined samples of |j0rgensen et al.ll2b0 6; Porras et al. 2007; 
iHarvev et al.ll200'7h . a small incidence of YSOs with protostellar 
disks misclassified as class I when seen edge-on could dramat- 
ically change the timescale of the embedded phase versus the 
protostellar disk phase. Taking just ~IQ% of stage II objects that 
present values of a typical of class I when seen edge-on, sim- 
ple algebra shows that a measured ratio of n(class I)/n(class II) 
= 26% would yield an intrinsic ratio of n(stage I)/n(stage II) 
= [n(class I)-n(edge-on disk)]/[n(class II)-i-n(edge-on disk)] - 
13%, half of the measured value. 

The effect of scattering is less strong in the silicate feature vs. 
the envelope mass (Fig.|4j;) than in the a panel. This is because 
silicates probe the entire column density, in contrast with ices 
which evaporate in the high density and high temperature inner 
part of the disk. Still, more massive envelopes produce larger 
values of silicate absorption than edge-on disks. This is not a 
column density effect but is due to the fainter overall continuum 
which in edge-on disks is increased by scattering. 
For the optical depth of the 3 fim H2O ice feature (Fig. |4}i) we 
also note that edge-on disks present values much lower than mas- 
sive envelopes despite higher column densities (see Fig.|2]i, again 
an effect of evaporation of ices combined with scattering of the 
surrounding continuum. A plot for the 15.2 fim feature of CO2 
is qualitatively similar to that for H2O with the exception that 
edge-on disks now reach values typical of 1.5 Mq envelopes. 



due to the smaller effect of scattering in the mid-IR. For the 
H2O ices there is no "traditional" boundary between class I and 
II. We adopted t - 0.5 because in this way the boundary crosses 
the Menv - 0.2 Mq models as it does for all the other indicators. 
In a similar fashion, below we chose the boundary for the CO2 
ice feature to be at t = 0. 1 . 

From the analysis of the results in Fig. H) it is evident that 
geometry has a different effect on the four indicators, so cross- 
correlating two evolutionary indicators could reveal different 
trends for the face-on and for the edge-on objects. This is shown 
in Fig.|5] where the silicate feature strength is plotted against a, 
Tboi, the optical depths of water ice and CO2 ice in panels a-d, 
while a versus Tboi is shown in panel e and the optical depths of 
the ices of water and CO2 are cross-correlated in panel f. The 
color coding and the shape coding are the same as in Fig|4] In 
each plot we added as small, black squares observations from the 
literature, in particular for sources studied by the c2d team with 
the Spitzer IRS. Values and references can be found in Table|2]in 
the online material. The models seem to reproduce fairly well the 
observational trends over the entire range of variables. Panels a-d 
show that indeed edge-on disks (open triangles) may not follow 
the same trends as the other models. If compared with embed- 
ded or face-on objects they systematically show deeper silicate 
absorption than would have been expected from the other indi- 
cators. As shown above this is due to the depth of the silicate ab- 
sorption feature being more sensitive to the total column density 
and less affected by scattering. This effect unfortunately is quite 
small and smaller than the current observational error on the in- 
dicators. The strong effect of scattering on the ice and silicate 
features makes it impossible to select the edge-on disks from the 
embedded YSOs from the different temperature structure along 
the line of sight. 

The Tbo\-a plane (Fig. |5^) is the one that is most populated by 
data. Most measurements come from the systematic work of 
Enoch et al. (in prep.) which makes use of 2MASS, Spitzer and 
Bolocam 1.1 mm data to calculate Tboi and a for all the YSOs 
in the Perseus, Serpens and Ophiucus star-forming regions stud- 
ied by c2d. The scatter in the observations is larger than in the 
models, which could be the result of foreground layers, incom- 
plete sampling of the SED, source confusion or the fact that in 
case of non-detections at 2.2 fim Spitzer bands were used to cal- 
culate a. Still, the general trend is reproduced reasonably well. 
Applying the effect of a foreground cloud (see arrows in Fig|5] 
and Sect. |3.4| i shows that while the sources in the bottom left 
corner can be still explained by these models, the sources that 
lie above the trend (l<Qf<3 and 100 K<rboi<600 K ) cannot and 
need further explanation. 

3.3. Sub-millimeter indicators 

An alternative to directly measure the envelope mass is through 
observation of optically thin tracers like the continuum at 
1.2 mm. In Fig.|6t, we show the predictions for the total 1.2 mm 
emission versus the mass of the envelope in our grid of models 
assuming an arbitrary distance of 140 pc. A correlation is seen 
when the envelope is more massive than 0. 1 Mq, but the depen- 
dence of the millimeter flux on the mass of the envelope becomes 
flat for stage II YSOs. This occurs because the disk starts to emit 
comparable amounts of flux. The single-dish flux can be cor- 
rected for the contribution from the disk by subtracting the flux 
measured in a much smaller region (Fig.|6]3), such as the one ob- 
tained by (sub)mm-interferometry. The ratio of single-dish over 
interferometer flux is perfectly correlated with envelope over 
disk mass (Fig.|6};), and offers a sensitive probe of evolutionary 
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stage independent of orientation. Using the visibility amplitude 
versus projected baseline length can further refine the distinc- 
tion between envelope and disk. With the increase in sensitivity 
of (sub)mm-interferometry currently reached at Plateau de Bure, 
SMA and CARMA, we expect these observations to be routinely 
done to check the true nature of class I objects. 



3.4. Discussion 

In our grid, we fixed some parameters that affect the density dis- 
tribution and so may affect the results. These were: the envelope 
centrifugal radius (fixed here at 300 AU), the disk and the en- 
velope outer radii (fixed at 200 AU and 10000 AU), the scale 
height distribution of the disk (fixed here to be 20% at 200 AU 
and decreasing with radius to the power of 2/7). Their effects 
are degenerate with those of changing the disk and the envelope 
masses, and for this reason we chose to keep them fixed. Here 
we explain how they change the results. 

Decreasing the centrifugal radius results in more peaked and 
more spherical envelopes, so for the same envelope mass we 
would get a larger column density. Varying the centrifugal ra- 
dius by a factor of 2 changes the column density by 30% in the 
same direction. Varying the disk and the envelope outer radii 
while keeping the same mass decreases the local volume density 
and thus decreases the column density along the line of sight. 
Doubling the outer radius reduces the column density by a factor 
of ~3, and vice-versa. Finally, varying the scale height distribu- 
tion may change the boundary inclination where, e.g., silicates 
turn from emission to absorption at a given mass. This boundary 
value was shown to range between 60° and 70° depending on 
the disk mass. An increase of 10% in the scale height reduces 
the boundary value by 5°; larger changes in the disk scale height 
are unlikely because of stability considerations. 

We can also calculate the effect of a foreground layer of a 
given column density, assumed to be dense enough to have ices 
throughout, by extinguishing the calculated SED and compar- 
ing with the original SED. As an example, having a foreground 
layer of 10^^ cm"^ (Ay ~ 5.5) would result in lowering a by 
0.21, dividing the silicate feature to continuum ratio by 1.46, 
and increasing the H2O and CO2 optical depths by 0.5 and 0. 1 
respectively. In the case of Tboi the effect of a foreground layer is 
neither linear nor logarithmic and depends on the details of the 
SED; on average, a Tboi of 1200 K would decrease by 350 K after 
being extinguished by a foreground layer of 10^^ cm"^ while the 
respective reductions for a Tboi of 600 K and 200 K are ~200 K 
and ~75 K. These reddening values relative to a 10^^ cm"^ fore- 
ground layer are represented with arrows in Fig.|5] 

Because of the rough parametrization of the density in the 
outflow cavity, we do not present any model results for face- 
on disks with / <25°. Test calculations show that widening of 
the outflow cavity to 25° has little impact on the observables 
for / >30°, but strongly affects them at /=0-20°. Consistent with 
previous work (Calvet et al. 1994; Whitney et al. 2003a) we find 
that Class I objects, seen below 15°, show flat spectra and sili- 
cates in emission, even for envelope masses up to 1 M©. 

We did not explore different (proto) stellar temperatures. For 
this reason, our models cannot be representative of high-mass 
YSOs with stellar temperatures higher than 5000 K, particularly 
in the case of the less embedded objects where the central radia- 
tion is n ot fully reprocessed (f or a better description of this effect 
see, e.g.. lWhitnev et alj|2004l) . 



4. Conclusions 

We used a Monte Carlo-based radiative transfer code to calcu- 
late the emission for a grid of models representing stage I and II 
young stellar objects: embedded YSOs and pre-main sequence 
stars with disks. The aim was to study the behaviour of classifi- 
cation parameters that are used to distinguish YSOs embedded 
in envelopes from those surrounded by disks only. We can sum- 
marize our results in these main points: 

- The classical indicators of YSO evolution, Q',rboi, and the sil- 
icate emission/absorption can describe very well the disap- 
pearance of the protostellar envelope for those systems seen 
at intermediate inclinations (approximately from 25° to 65°), 
accounting for more than half of the sources. Interestingly 
the transition between class I and class II measured with a, 
Tboi, or with the silicate feature happens for the same line 
of sight column density of ~ 1 ■ 10^^ cm"^ (corresponding 
to Menv~ 0.2 Mq), which means that these methods are well 
calibrated relative to each other. 

- Each of these indicators confuses pre-main sequence stars 
with edge-on disks with heavily embedded YSOs if the col- 
umn density exceeds 1 ■ 10^^ cm"^; this column density is 
reached for inclinations greater than 60°-70°, depending on 
the mass of the disk. 

- Even a pre-main sequence star surrounded by disk of only 
5 ■ 10""* Mq is classified as a flat-spectrum source for incli- 
nations greater than 70°. This means that 34% of the stage 
II sources are misclassified and could account for the entire 
population of flat-spectrum sources observed in nearby star- 
forming regions. If 10% of pre-main sequence disk objects 
are misclassified as class I, the derived time-scale for the em- 
bedded phase is overestimated by a factor of two, when using 
the ratio of class I/II blindly. 

- CO2 and H2O ice optical depth correlate well with a and 
Tboi, providing an alternative way to classify YSOs. For 
edge-on sources, the strong effect of scattering saturates the 
continuum around the features, thus affecting the feature 
strength and hiding the effect of temperature structure along 
the line of sight. 

- A combination of single dish and interferometer millimeter 
continuum observations is able to trace accurately, and with 
the minimum number of observations, the amount of remain- 
ing envelope in emission and thus the true evolutionary stage 
of the YSO. 
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Fig. 1. Density and temperature structures in two of our models. In the two top panels we show the density structure (logarithmic 
color-scale) and the dust temperature structure (black contours) of a model with a small envelope (Menv=0.074 Mq). The right and 
left panels show two different spatial scales. The two lower panels are similar to the upper two but for a model with a massive 
envelope (Menv=l-4 M©). From the comparison of the two models one can see that massive envelopes are also colder than smaller 
ones. 
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Fig. 2. Column density as a function of viewing angle for two of our models, one with Menv=0.074 M© (in the three top panels) and 
the other with Menv=l-4 Mq (in the three bottom panels). In all the panels, envelope material is represented with dotted Une while 
disk material is shown with solid lines. In panels a) and d) we show the total column density, while in b) and e) we show only the 
column density from cells whose temperature is lower than 90 K and in c) and f) only the cells with dust temperature warmer than 
90 K. Panels b) and e) thus represent the column density of material with ices, illustrating that for inclinations smaller than 65°-75° 
the envelope is the main contributor for ices along the line of sight while the disk dominates only at larger inclination. 
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Fig. 3. Normalized spectral energy distribution for a series of models with M(ji,sk=5 ■ 10"^ Mq and Lstar=6 L©. The six panels 
correspond to six different inclinations, and the 7 SEDs in each panel are for increasing Menv=0.01, 0.027, 0.074, 0.2, 0.55, 1.5 and 
4 Mq from top to bottom. The 10 fim silicate feature is in emission in the face-on objects with Menv< 0.07 Mq while it turns in 
absorption for inclinations greater than ~ 60° (see Sect. l3.ll l. In the same objects a, the spectral index between 2.2 and 24jt/m, turns 
from negative to positive, showing that both indicators fail to classify edge-on disks. 
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^envelope [^o] M^„^,l„p^ [MJ 

Fig. 4. Traditional class indicators versus envelope mass. Colors indicate envelope mass. If the system is seen face-on (filled circles), 
the envelope mass is well captured by a, Tboi, and the depth of ices and silicates. This is not true for edge-on systems (empty 
triangles), and in particular for edge-on stage II YSOs which are indistinguishable from stage I objects. To visualize better the edge- 
on from the face-on systems in these plots, edge-on systems are offset by increasing their envelope mass by 15%. The horizontal 
dotted line shows the traditional separation between class I and class II, while the vertical solid line denotes envelope masses of 
0.1 M0 separating stage I and stage II YSOs. Sectors with a discrepancy in classification between class and stage are hatched. 
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Fig. 5. Cross-correlation of various indicators. In panels a to d we compare the silicate feature depth with a,rboi, and the optical depth 
of water-ice and CO2 ice, while in panel e a is plotted versus Tboi and in panel f the two ices optical depths are cross correlated. In 
each panel observations from the literature are plotted as black squares. The boundary between class I and class II objects is shown 
with a vertical solid line and a horizontal dotted line. Face-on (filled circles) and some edge-on systems (empty triangles) form a 
well defined sequence that allow to relate the different evolutionary indicators. Color coding is the same as in Fig. |4] and indicates 
envelope mass. Besides the correlations in panels a to d, a slightly different trend seems to be found for edge-on disks, which present 
the heaviest geometrical effect. The effect on the evolutionary indicators due to a 10^^ cm"^ foreground layer are represented with 
arrows in each panel. Since the effect of the foreground on Tboi is not uniform we show the reddening vectors in 4 different ranges 
of Tboi (panels b and e) . 
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Fig. 6. Millimeter flux versus envelope mass. In panel a the envelope mass is seen to be correlated to the total flux at 1.2 mm 
(calculated assuming a distance of 140 pc). The envelope mass can be measured through this optically thin emission, but this can 
hardly be used to disentangle the edge-on disks from the low-mass envelopes, because in the transition between stage I and stage 
II, the masses of the disk and of the envelope become comparable, as are the two emissions (as shown in the right part of the plot). 
In panel b, we correct for the flux coming from the inner 300 AU (~2" at 140 pc) (obtainable with mm-interferometry) showing a 
much tighter correlation. In panel c we substitute the envelope mass with the ratio between disk and total mass and show that this 
correlates perfectly with the ratios between mm-interferometry and single-dish fluxes. 



Table 1. Observables for the models presented in Fig. [3] 
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= 5 ■ 10-' 


A^(H2) in disk [cm^^] 













3.61(16) 


5.56(20) 


1.07(22) 


1.76(23) 


2.75(25) 




= 0.010 


A'(H2) in envelope [cm" 


-2-1 


3.17(19) 


1.78(20) 
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Tm [K] 
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A'(H2) in envelope [cm" 
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a 
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1.72(21) 


2.09(21) 


2.49(21) 


2.22(21) 


2.68(21) 


4.09(21) 






a 






-0.65 


-0.66 


-0.68 


-0.71 


-0.55 


-0.20 


-0.48 






Tm [K] 






986 


981 


988 


984 


729 


251 


233 




= 0.20 


A'(H2) in envelope [cm" 


-2-1 


6.37(20) 


3.57(21) 


4.68(21) 


5.66(21) 


7.19(21) 


6.17(21) 


7.48(21) 


1.15(22) 






a 






-0.59 


-0.58 


-0.59 


-0.60 


-0.43 


-0.17 


-0.37 






Tboi [K] 






804 


767 


749 


720 


504 


172 


159 




= 0.55 


A'(H2) in envelope [cm" 


-2-1 


1.73(21) 


9.72(21) 


1.27(22) 


1.54(22) 


1.87(22) 


1.77(22) 


2.09(22) 


3.23(22) 






a 






-0.43 


-0.38 


-0.34 


-0.31 


-0.11 


0.026 


-0.082 






Tboi [K] 






486 


414 


374 


332 


224 


89 


82 




= 1.5 


A'(H2) in envelope [cm" 


-2-1 


4.71(21) 


2.64(22) 


3.46(22) 


4.17(22) 


5.01(22) 


7.02(22) 


5.81(22) 


8.50(22) 






a 






-0.038 


0.14 


0.26 


0.38 


0.59 


0.62 


0.64 






Tboi [K] 






188 


139 


117 


99 


75 


51 


48 




= 4.0 


A'(H2) in envelope [cm" 


-2-1 


1.28(22) 


7.18(22) 


9.40(22) 


1.13(23) 


1.35(23) 


1.74(23) 


1.74(23) 


2.38(23) 






a 






0.85 


1.3 


1.5 


1.6 


1.8 


1.8 


1.9 






Tboi [K] 






63 


51 


46 


42 


40 


38 


37 



Crapsi et al.: The true nature of class I YSO, Online Material p 1 



Online Material 



Crapsi et al.: The true nature of class I YSO, Online Material p 2 
Table 2. Observed indicators for various low-mass YSOs from the literature for sources in the c2d sample. 



Source 


^9.7 / CQiifinuum 


a(2.2-24fim) 




tHjO 


rCOj 


Ref. 


IRAS 04169+2702 


0.51 


1.29 


170 






1,3 


IRAS 04181+2654A 


0.68 


0.31 


346 






1,3 


IRAS 04239+2436 


0.52 


1.10 


236 






1,3 


Haro 6-1 OA 


0.41 


0.71 








1 


Haro 6-1 OB 


0.64 


1.56 








1 


IRAS 04287+1801 


0.29 


1.78 


97 






1,3 


LkCa 15 


3.18 


-0.90 








1 


IRAS 04381+2540 


0.39 


1.31 


139 






1,3 


Hen 3-600 A 


1.39 


-0.62 








1 


IRAS 14050-4109 


2.85 


-0.68 








1 


HD 163296 


2.33 


-0.58 








1 


HD 179218 


1.56 


-0.13 








1 


WW Vul 


2.65 


-0.57 








1 


HD 184761 


2.45 


-1.38 








1 


L1448 IRSl 


0.72 


0.34 




0.2 


0.064 


2 


IRAS 03235+3004 


0.17 


1.44 




3.7 


1.50 


2 


IRAS 03245+3002 


0.047 


2.70 






0.24 


2 


L1455 SMMl 


0.19 


2.41 




3.0 


0.90 


2 


RNO 15 


1.42 


-0.21 




0.22 


0.11 


2 


L1455 IRS3 


0.93 


0.98 






0.12 


2 


IRAS 03254+3050 


0.28 


0.90 




2.0 


0.39 


2 


IRAS 03271+3013 


0.39 


2.06 




1.9 


0.63 


2 


IRAS 03301+3111 


0.67 


0.51 




0.17 


0.028 


2 


Bl-a 


0.033 


1.87 






1.15 


2 


Bl-c 


0.005 


2.66 






1.95 


2 


Bl-b 


0.14 


0.68 






1.65 


2 


IRAS 03439+3233 


0.72 


0.51 




0.45 


0.10 


2 


IRAS 03445+3242 


0.36 


0.78 




0.80 


0.34 


2 


L1489 IRS 


0.36 


1.10 


238 


2.5 


0.65 


2,3 


IRAS 04108+2803B 


0.77 


0.90 


205 


1.4 


0.32 


2,3 


IRAS 04239+2436 


0.41 


0.79 


236 


1.2 


0.35 


2,3 


DG Tau B 


0.19 


1.16 




1.1 


0.24 


2 


HH46 IRS 


0.15 


1.70 




3.4 


0.85 


2 


IRAS 12553-7651 


0.31 


0.76 


99 




0.21 


2,4 


IRAS 13546-3941 


1.12 


-0.06 






0.40 


2 


IRAS 15398-3359 


0.048 


1.22 


61 




1.45 


2,4 


Elias 29 


0.27 


0.53 


409 


1.7 


0.16 


2,3 


CRBR 2422.8-3423 


0.23 


1.36 




2.1 


0.53 


2 


RNO 91 


0.55 


0.03 


715 


1.9 


0.45 


2,3 


IRAS 17081-2721 


1.37 


0.55 




0.55 


0.14 


2 


B59 YS05 


0.12 


2.26 






1.02 


2 


2MASSJ171 12317 


0.027 


2.48 






1.89 


2 


EC 74 


0.69 


-0.25 




0.41 


0.11 


2 


EC 82 


3.06 


0.38 






0.09 


2 


SVS 4-5 


0.11 


1.26 




2.8 


0.74 


2 


EC 90 


1.13 


-0.09 






0.21 


2 


EC 92 


1.31 


0.91 




0.68 


0.35 


2 


CK4 


1.45 


-0.25 






0.08 


2 


R CrA IRS 5 


0.43 


0.98 


736 


2.8 


0.69 


2,4 


HHIOO IRS 


0.52 


0.80 


373 


1.5 


0.24 


2,4 


CrA IRS 7 A 


0.091 


2.23 




4.1 


0.83 


2 


CrA IRS 7 B 


0.041 


1.63 




4.1 


1.10 


2 


CrA IRAS32 


0.093 


2.15 


148 




0.81 


2,4 


L1014IRS 


0.18 


1.28 




3.9 


1.02 


2 


IRAS 23238+7401 


0.20 


0.95 




4.0 


1.43 


2 


WL 12 


0.19 


1.75 


325 




0.19 


5,3 


WL6 


0.29 


1.28 


369 




0.40 


5,3 


IRS 63 


0.43 


0.57 


529 




0.31 


5,3 


IRAS 04302+2247 


0.25 


0.47 


202 




0.036 


5,3 


GSS 30 IRSl 


0.21 


1.30 


273 




0.13 


5,3 


L1535 IRS 


0.16 


1.26 


200 




0.25 


5,3 


2MASS 16282 


1.13 


-0.61 






0.00 
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